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^ The self-consistent harmonic oscillator model including the three- 

dimensional cranking term is extended to describe collective excitations in 
the random phase approximation. It is found that quadrupole collective 
excitations associated with wobbling motion in rotating nuclei lead to the 
appearance of two- or three-dimensional rotation. 
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Nuclear states grouped into AI = 2 sequences are remarkably well described 
in terms of the principal axis rotation which is the basis of the cranking 
model |1[. The principal axis cranking rotation was intuitively justified by a 
classical rigid body rotation which is favorable for a uniform rotation around 
the long or the short principal axis. However, if the classical body is not rigid, 
it may also uniformly rotate about an axis which does not coincide with the 
principal axes of the density distribution. In fact, more than a century ago 
Riemann |2| pointed out that such a situation could occur in the ellipsoidal 
self-gravitating fluid. The physical mechanism behind it is the vortex motion 
in a macroscopic system. 

Numerous experimental observations implying AI = 1 sequences in near 
spherical nuclei H raised the question about a new type of rotation and 
its physical nature in a finite quantum many body system. For particular 
nuclei the proton and neutron spin vectors may be oriented along different 
axes. Consequently, the total angular momentum could lie at an angle dif- 
ferent from one of the principal axes, i.e. at a tilted angle. It seems shell 
effects are the driving forces in this case which can be described within the 
Tilted Axis Cranking Model @] which assumes a two-dimensional rotation. 
Other intriguing examples are the rotational bands observed in Hf, W and 
Os (Z = 72 — 76) nuclei with A ~ 180 || |[ [/J. These bands are charac- 
terized by high K-values, where K is the angular momentum projection on 
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the body-fixed symmetry axis M. At finite rotational frequency a compe- 
tition could occur between states with large collective angular momentum 
oriented perpendicular to the symmetry axis and high-K states. Such com- 
petition gives rise to the new backbending phenomenon which is expected to 
be caused by the crossing of the ground and tilted rotational bands f§] . The 
one-dimensional cranking model based on the signature concept (the invari- 
ance of the wave function under rotation by tt around the rotational axis) 
fails in this situation, since signature is not conserved for tilted rotations. 
Furthermore, according to the standard Alaga rules fl], ||, the quadrupole 
B(E2) transitions to the ground band are forbidden due to the K selection 
rules (AK > 2). However, it is observed that some K-isomer states do decay 
via quadrupole transitions to the ground state. It was suggested that either 
fluctuations in the orientation of the angular momentum or shape fluctua- 
tions |6j or tilting degrees of freedom in the rotation axis H are responsible 
for this phenomenon. It was also observed that the AI = 2 sequence of the 
rotational states has been superseded by AI = 1 transitions in the yrast 
rotational band at high rotational frequency. In the present paper we con- 
siderably extend a previous study and demonstrate that tilted rotations 
naturally occur beyond a critical rotational frequency as an instability of 
collective vibrations caused by the wobbling motion. 

At high spin, the dynamical fluctuations of the shape and angular momen- 
tum vector can be described by the cranking random phase approximation 
(CRPA) as formulated in [|H], |TT| for a one-dimensional rotation. This ap- 
proach is the main theoretical tool used for the analysis of collective exci- 



tations in rotating nuclei (see [12|, [T3|, [14| and references therein). To un- 



derstand the main features of the formation of tilted bands, we assume the 
angular momentum and shape fluctuations to be dominant. For compari- 
son with experimental data pairing vibrations [0, |15| should be taken into 
account, but this degree of freedom does not change our main conclusions; 
it involves tedious calculations to be postponed for future publication. The 
CRPA and many RPA calculations for non-rotating nuclei suffer from the 
inconsistency between the basis generated by the mean field and the residual 
two-body interaction. To facilitate analytical and numerical results and to 
avoid spurious solutions due to such inconsistencies, we base our analysis 
upon the three-dimensional harmonic oscillator model. 
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We start with the single-particle Routhian 

- - N 1 m - - 

H n = H -n.L = E(^- 2 + + "ytfi + u\z]) - n • /,) (1) 

where the rotational vector Q of the cranking term has the components 
(Q x ,Q y ,Q z ) = f2(sin#cos0, sin # sin 0, cos 6). The Hamiltonian, Eq. ([!]), rep- 
resents the simplest mean field Hamiltonian which reproduces quite well the 
main properties of rotating nuclei || |TjJ . The spectrum of this Hamiltonian, 
i.e. the eigenmodes, is obtained by solving a third order polynomial in E 2 . 
The normal mode operators a\ and k = 1,2,3 ([afc,aj] = S^i) are lin- 
ear transformations in the coordinates and momenta pi, (i = x,y,z) with 
complex coefficients. In this way we obtain = Ej=i Ei=i,2,3 Ei(rii + 1/2) j 
with rii = ajaj. The normal modes are filled from the bottom which gives 
the ground state energy in the rotating frame E to t = E\ Ei +E 2 E2 +-^3 S3 
where =< Ej T ( n k + l/2)j > and the configuration is determined by 

Ei < E 2 < E 3 - 

The minimization of the Hamiltonian Eq.(0) with respect to the three fre- 
quencies u>i, subject to the volume conservation condition uj x uj y uj z = Uq, 
yields the self-consistency relation 

^ 2 >=^ 2 >=^V>. (2) 

These energy minima and their corresponding values for Ui depend on the 
rotational vector Q. For given Q we search for tilted solutions by seeking the 
minimum in the 9 — 0-plane ||. Three major results are reported here: 

1. starting from oblate, prolate or tri-axial shapes, for sufficiently high 
values of Q a local minimum, which is characterized by an oblate shape 
with the symmetry axis coinciding with the rotational axis, (9 = 90°, 
4> = 0°) is always obtained; 

2. increasing further the value of Q this local minimum shifts away from 
9 = 90° to 9 < 90° and possibly <fi > 0° thus giving rise to a two- or 
three-dimensional rotation; 

3. only when starting from tri-axial, near oblate, nuclei is this specific 
local minimum a global and stable minimum. For prolate nuclei, this 
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local minimum is attained only for values of fi so large that other (pos- 
sibly unphysical) minima with different configurations occur at a con- 
siderably lower energy, see Fig.l. When starting from the outset with 
an oblate nucleus, the global minimum is obtained for the rotational 
axis perpendicular to the symmetry axis; the specific local minimum 
referred to above is then the second lowest minimum with the same 
configuration. 

We first focus our attention to tri-axial near oblate nuclei. Here we encounter 
two critical values. At fig) there is a shape transition from tri-axial to 
axial (oblate) symmetry. This axial shape remains unchanged for a range 
^g"* — ^ — ^cr" 1 where fig) signifies the onset of a tilted rotation where 
the rotational axis no longer coincides with a principal axis. Note that for 
< fi < fig) the rotational axis is aligned to a principal axis which becomes 
the symmetry axis for fi > fig). At this point there is a confluence in 
the energy contours E(u y ,u z ) of two minima into one at u± = io y = uj z . 
In contrast, at the larger value fi = fig) the minimum at 9 = 90°, = 0° 
bifurcates into two minima as illustrated in Fig. 2. The two minima occurring 
at 9 = 75° and 9 = 105° (0 = 0°) appear at first glance physically equivalent. 
However, they are not as can be seen when fi is increased further in which 
case they move into positions with different values of (and 9) and only one 
is the global minimum. In fact, the directions of fi play a different role in 
the different octants [F7]. We interpret the rotation for fig) < fi < fig) as a 



K-isomer. For the tilted rotation the symmetric shape no longer prevails. 
While these findings are obtained from the numerical minimization procedure 
we now further substantiate our analysis by an analytical procedure. Using 
the fact that the first critical rotational frequency is associated with a one- 
dimensional rotation, we exploit in Eq.(|2]) that at the transition point uo y = 
u z = uj± and obtain a third order equation 

o u 1 r — 1 

it h = 3 

2 2r + l V ; 

where u = Q/uj± and r = X^/X^- From the discriminant of Eq.(^|) we 
obtain the critical value r cr = [\fTI + y / 2)(v / 27 — V2). It was shown in 



[ |18|j that a prolate system becomes eventually oblate with the rotational axis 
coinciding with the symmetry axis, if the initial deformation obeys r < r cr . 
However, according to our results, this case is not a global energy minimum, 
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in fact the global minimum is unphysical. In general, for r < r cr we obtain 
three solutions 

u = y| cos(^t^) n = 0,1,2 (4) 

21 + r < 5) 

These values of the rotational frequency correspond to three bifurcation 
points. Below we demonstrate that one of these solutions is the critical 
point 0^ where the lowest vibrational frequency tends to zero. 
The variations in the one-body potential around the equilibrium deformation 
determine the effective quadrupole-quadrupole interaction ([HJ. The total 
Hamiltonian can be presented as 

H = Hq — ~ £ QlQ, = H-n-L (6) 

fi=— 2 

While the mean field Hamiltonian Eq.(|l]) breaks the rotational symmetry, 
the Hamiltonian H in the total Hamiltonian Eq.(^) fulfills the commutation 



rules 



H, Li 



0, i = x,y,z. Note that the self-consistent condition Eq.(g) 
fixes the quadrupole strength constant k in the RPA calculations. We solve 
the RPA equation of motion for the general coordinates X x and momenta V x 
(see for details HI) 



[H,X x ] = -ilj x V x , [H,V x ]=iwxX x , [Xx,V x >} = i6 XjX > . (7) 

Here uo x is the RPA eigenfrequency in the rotating frame and the associated 
phonon 0\ = {X x —iP x )/\/2. In contrast to the CRPA approach, the phonon 
is in the present model a superposition of different signature phonons. The 
degree of the mixture depends on the tilted angle: the signature and |K| 
are good quantum numbers, respectively, for rotations perpendicular and 
parallel to the symmetry axis. The non-zero solutions appear in pairs ±hui x , 
we choose solutions with positive norm. 

Let us first focus on the RPA solution which leads to the transition from the 
non-axial to the oblate regime of rotation. Since at this transition point the 
minimal solution corresponds to a one-dimensional rotation (Q = (0, 0, 0)) 
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around the symmetry axis, the angular momentum becomes a good quan- 
tum number. The RPA states can be characterized by the projection of the 
angular momentum because [L x ,Oj] = XO\. Consequently, we obtain 



From Eq.(^|) it follows that at the rotational frequency Q cr = oj\/X one of 
the RPA frequency should be equal zero. The solution of the RPA equation 
for the positive signature quadrupole phonons with the largest projection 
A = — 2 gives 



This mode corresponds to the quadrupole de-excitation which leads to the 
state with two units angular momentum less then the vacuum state (K- 
isomer state). From Eq.® we obtain the first critical value of the rotational 
frequency at which the transition from the non-collective rotation to the 
non-axial collective rotation takes place 



Therefore, the positive signature quadrupole de-excitation leads from the 
K-isomer state to the yrast state with non-axial quadrupole shape. 
Quantum excitations describing the wobbling motion correspond to the 
negative signature quadrupole phonons [|T0| , [TTj |20| . These excitations are 
connected to the yrast line by means of the quadrupole transitions which 
carry one unit of the angular momentum. Similar to the positive signature 
phonons, the negative signature phonons can have zero excitation energy at 
particular rotational frequency. Solution of the RPA equation for the nega- 
tive signature phonons with A = —1 for the oblate rotation regime leads to 
the result 



[H, 0{] = [H- QL X , 0{] = (u)\ - \n)0[ = u x O 



(8) 






(10) 




(11) 



cosifj = 3\/3 



r — 1 



(12) 




,2 )3/2 l + r 
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The condition lox=^i = yields the second critical frequency 



«« = \/S (cos |_V38i4). (13) 

Beyond this rotational frequency the mean field solution corresponds to the 
stable tilted rotation. We recall that the expressions given in Eqs. fllOD and 
( |13D coincide exactly with the transitional points found in the numerical 
minimization procedure. Note that, in contrast to the familiar phase transi- 
tion from spherical to deformed shape owing to the variation of a strength 
parameter [|lj , it is here the rotation that leads to a phase transition from tri 
axial to oblate shape and then to the tilted rotation associated again with a 
non-axial shape for the same strength parameter of the residual quadrupole- 
quadrupole interaction. 

According to our analysis, only nuclei near to oblate shape, when non- 
rotational, exhibit at a certain rotational frequency tilted rotation which, 
for lesser rotational speed, leads to the high K-isomer states. These states 
decay through non-axial shapes to the ground states. Nuclei which are oblate 
when non-rotational also can have this behavior which corresponds to the 
second (local) minimum discussed above. However, they also have states of 
good signature with an even lower energy for the same rotational speed and 
the same configuration. In contrast, for prolate nuclei the local minimum 
with these characteristics occurs only at a value Q so large that the global 
minimum occurs for an absurdly unphysical configuration; for this reason it 
is doubtful to attach particular physical significance to the local minimum. 
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Figure Captions 

Fig.l Energy contours in the uj x — u) y — plane for a rotational speed where 
the onset of an oblate shape has just taken place and corresponds to the 
local minimum shown in the top. However, for the nucleus considered (N = 
44, tri-axial near prolate) a much deeper minimum has developed as an 
instability as seen in the top right corner of the bottom picture. The section 
of the top illustration is indicated in the bottom part. 

Fig. 2 Energy minimum in the — #-plane just before (top) and just after 
(bottom) the onset of tilted rotation. Note the local maximum developing at 
the centre. The example chosen is N = 36, a tri-axial near oblate nucleus. 
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